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ABSTRACT 

We investigate the acceleration and simultaneous radiative losses of electrons in the 
vicinity of relativistic shocks. Particles undergo pitch angle diffusion, gaining energy as 
they cross the shock by the Fermi mechanism and also emitting synchrotron radiation 
in the ambient magnetic field. A semi-analytic approach is developed which allows us 
to consider the behaviour of the shape of the spectral cut-off and the variation of that 
cut-off with the particle pitch angle. The implications for the synchrotron emission of 
relativistic jets, such as those in gamma ray burst sources and blazars, are discussed. 
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1 INTRODUCTION 

The role of radiative losses in determining the spectra from 
non-thermal sources has been well und erstood in the non - 
rela tivistic shock limit since the wor k of I Webb" et all (|l984l ) 
and lHeavens fc Meisenheimerl (|l987l ). Their results were in 
broad agreement with the natural expectation that there 
would be a cut-off in the spectrum, at the shock, and at a 
momentum where acceleration and loss timescales are equal, 
with the shape of this cut-off depending critically on the 
momentum dependence of the particle scattering. Subse- 
quently, as the particles are advected downstream, and are 
no longer efficiently accelerated by the shock, the spectra 
steepens at momenta where the particles have had sufficient 
time to cool. At a strong, nonrelativistic shock the differen- 
tial number density of particles at energies where radiative 
cooling is unimportant is a power law with N(E) oc E~ a ' 5 
with a corresponding intensity of I v oc v~ ' 5 for the emitted 
synchrotron radiation. At higher momenta, where cooling 
becomes important, the spectrum steepens so that the ra- 
diation, beyond a break frequency Vb, is I v oc v^ 1 up to a 
critical frequency, v c corresponding to cut-off of the particle 
spectrum. The position of v\, depends on position away from 
the shock; decreasing downstream as the particles have more 
time to cool. The observed emission is therefore dependent 
on the spatial resolution with which the sourc e is observed 
as discussed in lHeavens fc Meisenheimerl (| 1987h . The results 
in the existing literature refer only to non-relativistic flows 
and are of great use in analysing the spectra from super- 
novae and the jets of some active galaxies (AGN). How- 
ever, a number of objects of astrophysical importance, such 
as AGN jets, microquasars and gamma-ray bursts, contain 
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flows which have bulk relativistic motion and the purpose 
of this paper is to examine the breaks, cut-offs and emission 
for such sources. 

While the first order Fermi process at relativistic shocks 
contains the same basic physics as in the nonrelativistic case, 
i.e. scattering leading to multiple shock crossings competing 
with a finite chance of escape downstream, the anisotropy 
of the particle distribution comp li cates the analysis c onsid - 
era bly dKirk fc Schneider] l|l987f ). iHeavens fc Drurvl (|l988l ) 



and iKirk et al.l ( 20001 )). The inclusion of self-consistent syn- 
chrotron losses will, as in the nonrelativistic limit, modify 
the spectrum at high momenta but we would also expect 
pitch angle effects to become apparent in the position of 
the cut-off and the emission itself. In order to motivate our 
treatment of this problem we first discuss the nonrelativistic 
shock limit in section 2, including the emission from a spa- 
tially integrated source. Section 3 then presents the analysis 
of synchrotron losses at relativistic shocks with particular 
emphasis on the shape of the momentum cut-off. We con- 
clude with a discussion in section 4. 



2 NONRELATIVISTIC SHOCKS 

The effect of synchrotron losses on the energetic particle dis- 
tribution in the pres ence of nonrelativis tic shocks is demon- 
strated rigorously in IWebb et al.l (Il984l). However a si mpler 
approach is described in IHeavens fc Meisenheimerl (|l987l ) 
provided synchrotron losses are not considered important 
at the injection energies. We will follow this approach here, 
although we shall introduce a slightly different definition of 
the cut-off momentum. 

In the presense of a magnetic field charged particles 
emit synchrotron radiation with an energy loss rate given 
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by 



dp 
~dt 



-a s B 2 p 2 = — Xp 2 



(1) 



where a s is a positive constant. The radiative loss timescale 
is therefore ti oss = l/(Ap). In the steady state, and in the 
presence of a nonrelativistic flow it, energetic particles obey a 
transport equation describing advection, diffusion, adiabatic 
compression and radiative losses, 



df 



— f K—— 

dz dz V dz 



1 dudj_ 

3 P dz dp 



i_d_ 

p 2 dp 



(Ap 4 /)=0. (2) 



In the presence of a nonrelativistic shock front where 
the upstream flow speed is it_ and that downstream is u+ 
the acceleration timescale is 

3 



^acc — 



(3) 



At momenta for which i acc <C t\ OBB the phase space den- 
sity will be a simple power law with / oc p~ 3 where 
s = 3u-/(u- — u+). 

2.1 Momentum Cut-off 

The spectrum will steepen at momentum p* where 
tacc(p*) = ti OSB (p*). In the case of momentum independent 
diffusion this gives 



u+ 



X-K- 



+ 



U+ 



(4) 



In the case of a relativistic shock this result no longer 
strictly holds since the acceleration timescale defined above 
is only valid for nonrelativistic flows. Nevertheless, we will 
use this definition of p* throughout the paper for the sake 
of comparison. 

However, we require a general definition of the cut-off 
momentum that can be applied in the relativistic limit. An 
obvious alternative is to define the momentum at which the 
local spectral index, 91n//<91np, becomes s + 1 but, as we 
shall see, it is necessary to perform a Laplace transform of 
the transport equation to proceed with this problem and 
it is more straightforward to define the cut-off in terms of 
spectral steepening of the Laplace transformed spectrum. 
In order to motivate such a definition we solve the nonrel- 
ativistic shock acceleration problem in the presence of syn- 
chrotron losses by first making the substitutions W = p 4 / 
and y = 1/p so that the transport equation, either upstream 
or downstream of the shock where adiabatic losses are zero, 
becomes 



dW d 



dW\ d .,„,. „ 
K ^J + ^ (W) =°- (5) 



dz dz 

Taking the Laplace transform with respect to y 

W(k,z)= / W(y,z)en>(-yk)dy (6) 



and using the fact that losses prevent any particles achieving 
infinite energy, i.e. W(0, z) = 0, the transformed transport 
equation is 



dW 



d dW 
dz \ dz 



+ XkW = 0. 



(7) 



in the case of a momentum independent diffusion coefficient. 



Since the distribution function must be bounded infinitely 
far upstream and downstream, the solution becomes 



W± = A±(k) exp I — z 



where we have introduced 



4A±k± 



(8) 



(9) 



The isotropic and anisotropic parts of the particle dis- 
tribution function must match up at the shock giving, 



/_(p,0) = /+(p,0) 



(10) 



df- . u- df- df+ u+ df+ 

dz 6 dp dz 3 dp 

Multiplying the isotropic boundary condition by p 4 , making 
the substitutions as above and taking the Laplace transform 
with respect to 1/p gives 

W-(k,z = 0) = W+(k,z = 0) (12) 

which in turn gives 

A-(k) = A+(k) = A(k). (13) 

The flux continuity condition (|11[) becomes 

A(k) = A k 3 ~ s exp (s\/l + u-k- (s- 3)\/l + uj + kj x 

+ ui+k) 

(14) 

In the absence of synchrotron losses, u± = 0, we have 
W±(k) oc fc 3_s which, upon inversion, gives f(p) oc p~ s as 
expected. We can therefore define a function, Q, by W = 
k 3 ~ s Q. Recalling that k is the Laplace transformed variable 
of inverse momentum we define the cut-off momentum, p C ut, 
to occur at the point where 



ainQ 



dink 



(15) 



k—p c 



As an illustrative example, consider a power law distribution 
with a sharp maximum momentum, f(p) oc p _s i/(p ma x — p) 
with H the Heaviside function. In this case we have 



W = y s ~ 4 H(y 



1 



Pmax 



W 



T(s -3)k 3 



(16) 



where F is the Gamma function. With Q oc exp(— fc/p max ) 
we then have p cut = p max as required physically in this sim- 
ple case. 

Returning to the solution of the shock problem we have 
from equation 1141 



= exp y—s\/ 1 + uj-k — (s — 3)y/T+ to+kj 



1 + y/l + uj-k) (1 + yjl+u+k 



Defining 

x = VW- 



(17) 
(18) 



2 s 2 n + n 2 s 2 



x (s - 3) 



6 2 s + 2Qs + 9fi 2 -8fi) 
(19) 
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le-04 0.01 1 100 _ 10000 

Figure 1. The cut-off momentum, p cu t as a function of f2 for 
fixed equilibrium momentum, p* = 1, and spectral index, s = 4. 



o.i 



0.01 



p* 




p 

r cut 


Sillzer Sum 






\l 1 






M=2 






M=6 







0.01 0.1 1 , 10 

log p 

Figure 2. Laplace inversion for s = 4 and Q = 1. Using just M = 
6 in the Salzer summation the inversion has already converged. 
The first approximation M = 1 is exactly the Laplace function 
Wk. We can see how fast the Salzer summation Post-Widdcr 
inversion converges as M = 2 is a very reasonably approximation 
to the actual solution. 
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Figure 3. Particle Distributions of Various Q 



n 


V* 


Pcut 


P 





1 


1.25 


2.25 


i 


1 


1 


2.25 


9 


1 


1.4 


2 


16 


1 


1.53 


1.8 


25 


1 


1.63 


1.75 


oo 


1 


2 


1.5 



Table 1. Parameters for fitting particle spectra 

to the distribution of the form 

exp (- (p/j5 cu tf ) (22) 

where (3 ~ 2. Table [1] shows how (3 varies with ft for a 
shock of natural spectral index s = 4, with f3 attaining its 
maximum value of 2.25, i.e. the cut-off is sharpest, when ft = 
1. When ft 3> 1 particles can diffuse in the upstream without 
losing any energy, allowing a greater spread in momentum 
above p cut . 



gives 



s 4 (2 + 2fi) - s 3 (5 + lift) + s 2 (5 + m - 2 X ) 



+ 



(s 2 (l - ft) +6fts 

s(33ft + x) ~ 36ft 
(s 2 (l - ft) +6fts - 9ft) 2 



9ft)' 



This is always greater than p* , 



(20) 



2.2 The Integrated Distribution Function and 
Synchrotron Spectra 

When the source cannot be fully resolved observationally, 
we must include the contribution from all particles within 
some distance z of the shock in calculating the spatially 
integrated emission. In the case of steady emission from a 
jet pointing towards us, or a completely unresolved source, 
z' is essentially the source size in the optically thin limit. For 



uj- s + (s — 3)ft 



(21) 



as can be seen from figure [T] The minimum value for p cut 
occurs for ft = 1, and is given by p cut = 2§^3 P*- 

The Laplace inversion (see Appendix for details) for s — 
4 and ft = 1 is shown in figure [2] Using just M — 6 in the 
Salzer summation the inversion has already converged. The 
first approximation M — 1 is exactly the Laplace function 
Wk. We can see how fast the Salzer summation Post-Widder 
inversion converges as M — 2 is a good approximation to the 
actual solution. 

Figure |3] shows how the particle distribution varies with 
ft; if ft 7^ 1 the cut-off is broader. While p* is independent 
of ft, pcut increases as the distribution broadens. 

The data in figure[3]can be fitted by an exponential tail 



logp 4 f 




Z '=5d z'=ia z'=i , \ z'=o 

0.01 0.1 1 log p 10 

Figure 4. The spatial variation of the particle distribution for 

n = i. z > = ^+ z . 
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Figure 5. The spatially integrated particle distribution for Q = 1 
from z' = to z' = 50. Note that as well as the particle cut-off 
there is a spectral break earlier were they spectrum softens from 
p~ 4 to p~ 5 . This is analogous to synchrotron ageing. 



simplicity we assume that the magnetic field downstream of 
the shock is constant although the model can be generalised 
for more complex cases. 

The integrated Laplace distribution function is 



/' 

Jo 



W dz 



A k + 



(1 + y/l+LU+k) 



11+ 

X 



uj+k 



exp 



Jl + U! + k u+ i 

2 

2 K+ 



(23) 



When z' is very small the result is k 3 ~ s with a cut-off at 
high k as expected. As z' tends to infinity at low k we have 
k 2 ~ s so the spectrum is steepened , with the same high k 
cut-off. For finite values of z' the spectrum starts as k s ~ s 
before turning into k 2 ~ s and finally cutting off. We shall see 
later that this result also holds in real momentum space. 

Figure [4] shows how the cut-off tends to lower momenta 
as we go further downstream. However what is most often 
observed a result of the integrated distribution is shown in 
figure [5] While the cut-off momentum is the same as at 
the shock, the distribution changes from an initial p~ 4 to 
a p~° spectrum at some critical momentum, ph, which de- 
pends on z' . Here we will consider only synchrotron emission 
from an ordered magnetic field (parallel to the flow). Let 

w = 47r ^^ c where v is the frequency, q is the charge on 
the electron, m e is the electron mass and B is the magnetic 
field strength. Then given a spatially integrated particle dis- 
tribution / o c p~ s g(p,n) the total pow er emitted per unit 
frequency is (|Rvbicki fc Lightmanlll98q ) 



Ptotiuj) OC \Jl - fl 2 




5)/2 



where F is the first synchrotron function 

f oo 

F(x) = x Ks(y) dy 

J X 



H F(x) dx (24) 



(25) 
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Figure 6. A telescope with very high resolution may be able 
to observe synchrotron radiation at the shock. The only spectral 
feature here is the cut-off at v C ut- Before u cu t the spectrum has 
shape I(u) ~ v~ x ^ 2 . For this result fi = 1. 
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Figure 7. Spatially integrated emission for £7 = 1. The black 
dashed curve shows the synchrotron spectrum of an unresolved 
object containing a strong shock. The only feature is the spectral 
cut-off at i^cut- Before i/ cu t the spectrum has shape I(u) ~ v^ 1 . 
The solid curve illustrates the synchrotron spectrum of a partially 
resolved source with emission from the shock (z' = 0) to some 
downstream distance (z' = 50). 



tion / downstream of a shock is assumed to be isotropic in 
which case g is independent of /i. 

Figure [5] shows the emission at the shock. The only fea- 
ture here is the cut-off hump, which is of course related to 
the particle cut-off, before which /„ ~ u~^ a ~ 3 ^ 2 . Figure [7] 
shows a more realistic plot, that of emission from an ex- 
tended region. If the region is infinite in size then the cut-off 
remains the sole feature but the spectrum before the cut-off 
is different I v ~ v~^ s ~ 2 ^ 2 . If the region has finite size then 
a second feature, the spectral break vt, appears. Before the 
break the spectrum goes as /„ ~ ^-( s - 3 )/ 2 w hile after it it 
Again this is related to the momentum 



is 



-C«-2)/2 



break ph we see in the particle distribution in figure [5] 



3 RELATIVISTIC SHOCK ACCELERATION 
WITH LOSSES 

In the case of a relativistic shock, the particle transport 
equation describing advection, pitch angle diffusion and 
losses becomes 



r (u + fi) 



'11 



dfx 



Dn 



2i 



1 d( P 4 f) 
p 2 dp 



(26) 



In the case of non-relativistic diffusive shock accelera- 



where fj, is the cosine of the pitch angle of the particle 
and the flow velocity is constant upstream and downstream 
of the shock. A = 2 % T ~ 2 Ub and g(fi,) = 1 — (i 2 for syn- 
chrotron losses in an ordered magnetic field, A = 4<ttc/b/3 
and g(fj.) — 1 for synchrotron losses in a tangled magnetic 
field, or A = 4<7T'7 ra d/3 and g(fj,) = 1 for inverse Compton 
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losses. Equation p6[) holds separately upstream and down- 
stream with the conditions that the distribution is isotropic 
infinity far downstream, there are no particles infinitely far 
upstream and the distribution is continuous at the shock. Al- 
though we will derive equations for general momentum inde- 
pendent pitch-angle diffusion and an arbitrary magnetic field 
alignment, the figures and results produced throughout the 
rest of this paper are for isotropic diffusion D^^ = D(l — fi 2 ) 
in an ordered (longitudinal) magnetic field with X/D — 0.1. 

Guided by the treatment of the nonrelativistic case we 
set W = p 4 f and y = 1/p so that 



r(w + /i)— =— 

az ofi 



o dW 



dW 

AffOx)-^. (27) 



Taking the Laplace Transform with respect to y and assum- 
ing W(0,ii,z) = 



TV j. \ 9W 9 



o dW 



Xg(n)kW. (28) 



With the spatial and pitch angle variables separable we look 
for solutions of the form 



(29) 



putting this back into the reduced transport equation we get 



ay 

r(u + fj,)-^-Qi = (VQi)Xi 
az 

where we have defined the differential operator T> via 



254 = #" 
on 



D( M )(1- M 2 ) 



Separating X and Q we get the usual 



- \g(fi)k$. 



i(u + fj.) 
Ai{k)z 



=> Xi{k, z) — exp 

and we have an equation for Q(k, fj,) 

VQi - K l {k)Q l {u + (i) = 0. 
Expanding out the differential operator we get 



t)(t 



(30) 

(31) 

(32) 
(33) 

(34) 

= 
(35) 



which has regular singularities at n = ±1 and so it should 
be possible to find solutions for Qi on [—1, 1] for all k G C. 



3.1 Determining the Eigenfunctions 

We know that along the real axis, k — x £ R, each Qi 
satisfies 



0_ 



D(,)d-, 2 )f 



(Ai(u + n) + xAg(n)) Qi = 0. 



(36) 



We define an inner product by: 

(C,0 = / (u + riCtdfi. (37) 



u_=.3; ^=.076 
u_=.5; u+=.129 ■ 
u_=.7; u^.189 



Figure 8. The zeroth order downstream eigenvalue for shock 
speeds it_ = .3, .5, .7. Along the x-axis we have plotted the log- 
arithm of fc_|_ while along the j/-axis we have the logarithm of 
— Ao(fc+). When k + = we have Aq(0) = 0. 



It can be shown that the Qi(x) are orthogonal and either 
real or purely imaginary, and the At (a;) are real and distinct. 
We can normalise the eigenfunctions such that 



(Qi, Qj) = 6i,j 
or considering them as real 

{Q i ,Q i ) = Si, i (l/2-i)/\l/2-i\=rn, 
Then we have (see Appendix for details) 

and 

^r 1 = -A ( 



(38) 



(39) 



(40) 



(41) 



We solve equation | 36l at x = using the Priifer transforma- 
tion as in lKirk et al.l l|2000t ). We then use equations 1401 and 
1411 to find Qi(x,/j,) and Ai(x) for x > using Runge-Kutta 
methods. 

Figures [S] [TO] and QT] show the zeroth downstream 
eigenvalues and eigenfunctions for shocks speeds of .3, .5 
and .7. This eigenfunction is the dominant component in the 
downstream distribution function at the shock of such mildly 
relativistic shocks, where we are close to isotropy. Further 
downstream, where the contribution of higher eigenfunctions 
are more strongly damped, so the anisotropy for some z > 
is essentially that of the zeroth eigenfunction. While Ao(fc+) 
is initially zero, note from figure [S] that it decreases linearly 
until a certain point which, as we will see later, is close to 
the cut-off momentum. This will play a major role in the 
integrated distribution function and emission. Figures |§J 
1101 and QT] show how anisotropy arises in the zeroth order 
eigenfunction which is isotropic for k — 0. Given that k is 
related to the momentum these figures show that, since this 
is the dominant eigenfunction, the anisotropy will increase 
with increasing energy. 
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Figure 9. The zeroth order downstream eigenfunction for shock 
speed U— = .3. Along the x-axis we have plotted the logarithm 
of fc+ while along the y-axis we have fj,+ . Up the 2-axis we have 
plotted Qo(^+i which also defines the grayscale. Note the 

anisotropy increases with fe+. 
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Figure 10. The zeroth order downstream eigenfunction for shock 
speed u— = .5. Along the x-axis we have plotted the logarithm 
of fc+ while along the y-axis we have fj,+. Up the 2-axis we have 
plotted <2o(&+> which also defines the grayscale. Note the 

anisotropy increases with fc_|_. 



3.2 Shock matching conditions 

Starting from 

W(k, fj,, z) = ai(k)Xi(k, z)Qi(k, fi) 

i 

^>(fc)exp^k 



(42) 

Qi(k,fj,)) (43) 



we note that upstream (z < 0) we have = for all i 
such that A, ^ and that downstream (z > 0) we have 
a* = for all i such that Aj > 0. The distribution function 
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Figure 11. The zeroth order downstream eigenfunction for shock 
speed U— = .7. Along the x-axis we have plotted the logarithm 
of fc+ while along the y-axis we have fj,+. Up the 2-axis we have 
plotted Qo(k+, fJ.+ ), which also defines the greyscale. Note the 
anisotropy increases with k+. 



is continuous at the shock, 

/-(j/-.M-,0) = f+(y+, n+,0) 



(44) 



with (y_,^_) related to (y+,fi+) by a Lorentz transforma- 
tion of velocity u ro i = (u_ — u+)/(l — U-U+), 

V~ = r ro iy + (l + Wrel/U-). (45) 

In terms of W the matching condition becomes 

it , ( 1 + u rclM _ fW- (y- ,^,0) = W+(y + ,n + ,0) (46) 

and we now need to express this in terms of W, the Laplace 
transform with respect to y. Taking k-y- = k+y+, multi- 
plying the matching condition for W by exp(— k+y+) and 
integrating over y+ gives 



rL(l + UreW-) 3 W-(k-, n-,0) 



W+(k+,v + ,0). (47) 



Guided by the discussion for the nonrelativistic case, we 
use the expansion 



ir d 



J2btkr +3 Qt(k±,V±) (48) 



so that the matching condition for the Laplace transformed 
spectrum at the shock reduces to 



r? e i(l + UydH-Y^bi (k-)Qi (k-,H-) = 

i 

j2t>i(k + )Qi(k + , f i+) 



(49) 



In order to solve for the particle spectrum, we multiply by 
(tt+ + A t +)Qj h (^+i M+) 3 ^ and integrate over yit+. Then for 
a fixed k- we have 



(u+ + ^+)Q+(fc + ,/i + )d/i+ = 0. 



Defining a matrix S with elements 



(50) 



(l + Urelfi-YQi (k-,(J,-)(u+ +/i+)Q+(fc + ,^ + )d^ + 

(51) 

we need to find the spectral index s, such that det S = 0. 
The Laplace inversion is then carried out numerically (see 
Appendix for details). As motivated by the nonrelativistic 
case, we define the cut-off to be the point at which 



djlnR) 
d(lnfc) 



(52) 



where R = J2 b iQi- Figures \T2\ UM and [14] plot 
rf(ln R)/d(\n k) at the shock against k as measured down- 
stream. The results are summarised in table [2] 

Figures [121 1131 and 1141 show how the cut-off momentum 
becomes increasing anisotropic as the shock speed increases. 
The distribution can be fitted approximately by 



exp 



P 



TrolPcut (AM 



(53) 



where /3 is typically 2. However it is difficult justify the use 
of the factor \JT\ci in general as our results are only for 
mildly relativistic shocks. This fit justifies our definition of 
p cut instead of using the equilibrium momentum p* . Figure 
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Figure 12. Plotted along the x-axis we have the logarithm of 
momentum fc+ while along the j/-axis we have d(ln i?-|_)/d(ln fc+) 
for «_ = .3 and = J]V b^(fc_|_)Q^~(fc+, Note the cut-off 
depends on /x+ 
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Pcut (M+ 
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Pcut (A*+ 


= 0) 


.64 
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Pcut (A*+ 
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Pcut 
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.741 
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Table 2. Summary of Cut-Off Momenta 
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Figure 13. Plotted along the x-axis we have the logarithm of 
momentum fc+ while along the j/-axis we have d(ln iJ+)/d(m fc+) 
for ii_ = .5 and R+ = J^. bf (k+)Q'£'(k+, fj,+ ). Note the cut-off 
depends on fi+ 



Figure 15. The downstream function evaluated at the shock for a 
shock speed of .3. Along the x-axis we have plotted the logarithm 
of momentum p+ while along the y-axis we have the logarithm of 
W = p 4 f. 



1171 illustrates this approximation for a .7c shock. j3 seems 
to be pitch angle dependent varying between 1.75 and 2.2, 
but typically 2. In fact for the .3c and .5c shock cases j3 
showed much less variation about 2. For the shock speeds 
we have chosen, with the Juttner-Synge equation of state, 
the spectral indices in the absence of losses are close to 4. 

Figures 1151 [161 and [17] illustrated a feature that was 
not present in the non-relativistic case. The pitch angle de- 
pendence of the cut-off momentum leads to a difference in 
the isotropy levels between particles above and below some 
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Figure 14. Plotted along the x-axis we have the logarithm of 
momentum k+ while along the y-axis we have d(ln i?+)/d(m fe+) 
for u- = .7 and R + = J2i b t ( fc +)Qi"( fc +i /•*+)• Notc tnc cut-off 
depends on /j+ 



critical momentum p a niso- Indeed there is a clear pattern of 
greater levels of anisotropy at high energies as the shock 
speed increases, despite the fact that the results presented 
here are only for mildly relativistic shocks. 

3.3 The Spatially Integrated Distribution 

While the method we follow in this paper finds the upstream 
particle distribution directly, it is easy to find the down- 

i , , — , — , — , — I 
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Figure 16. The downstream function evaluated at the shock for a 
shock speed of .5. Along the x-axis we have plotted the logarithm 
of momentum p+ while along the j/-axis we have the logarithm of 
W = p 4 f. 
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Figure 17. The downstream function evaluated at the shock for a 
shock speed of .7. Along the x-axis we have plotted the logarithm 
of momentum p_|_ while along the jy-axis we have the logarithm 
of W = p 4 f. The lines are data while the points are the best fit 
described in the text. 
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Figure 18. The downstream function integrated for a shock 
speed of .3 between z' = and z' = 100 where z' = Dz/T-l.. 
Along the x-axis we have plotted the logarithm of momentum p_|_ 
while along the y-axis we have the logarithm of W = p 4 f. 
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Figure 19. The downstream function integrated for a shock 
speed of .5 between z' = and z' = 100 where z' = Dz/T-l.. 
Along the x-axis we have plotted the logarithm of momentum p_|_ 
while along the j/-axis we have the logarithm of W = p 4 f. 
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Figure 20. The downstream function integrated for a shock 
speed of .7 between z' = and z' = 100 where z' = Dz/F+. 
Along the x-axis we have plotted the logarithm of momentum p_L 
while along the j/-axis we have the logarithm of W = p 4 /. 



stream distribution by using the matching condition, as dis- 
cussed in the previous section. The downstream distribution 
is. in many respects, more important physically as it will 
be responsible form most of the spatial integrated emission. 
As it can be difficult to spatially resolve observational data 
from non-thermal emitters, we must consider the emission 
from an extended region of space. Our eigenfunction expan- 
sion allows us to do this quite easily. The spatially averaged 
distribution from a downstream region [zo,zi] in terms of 
Laplace variables is 

/ ( A-t(k+)zi 
%o,*i]( fe +>/ i +) = Z^ a *( k ) ( ex P ( p 

exp (Ai(k + ) z ^y + QU, + , k+ ), (54) 



In the case of a source which is completely spatially unre- 
solved this reduces to 

fl [0iOol (k+ ,/,+) = -r + ^ <H(k) Qt }^' \ +) . (55) 

Of course the optical depth of the emitting region will also 
have an effect on the spectrum of unresolved sources by re- 
ducing Z\. 



Using the same numerical Laplace inversion as in the 
non-relativistic case we have calculated the distribution 
functions and synchrotron emission. Figures 1181 1191 and [20] 
show the spatially integrated distribution functions for a fi- 
nite emission region. Now there are there features: a momen- 
tum break, pb, due to spatial effect; an anisotropic break, 
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Figure 21. Synchrotron emission from the particle distribution 
shown in 1181 measured in the downstream medium. In plotting 
our fi = ±1 we used fi = ±.9999 as there is no emission from an 
ordered field along fj, = ±1. 
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Figure 22. Synchrotron emission from the particle distribution 
shown in 1191 measured in the downstream medium. In plotting 
our fi = ±1 we used fi = ±.9999 as there is no emission from an 
ordered field along fj, = ±1. 
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anisotropy emerging resulting in synchrotron emission which 
is also anisotropic. This could be extremely important in the 
modelling of th e inverse Compton hum p in 7-rays observed 
111 TeV Blazars {Ah aroman et alj2006h . As a second implica- 
tion of the particle anisotropy, in the presence of losses, the 
idealised situation, of a two sided strongly polarised iden- 
tical jet system can be considered. Each jet contains only 
forward external shocks, and the jet which is directed to- 
wards the observer is inclined at an angle 9 — cos _1 (— fi) 
to the line of sight (magnetic field direction same as that 
of shock). Then we will observe the emission from particles 
in the jet directed towards us which have pitch angle fi and 
from particles in the jet directed away from us which have 
pitch angle — fx. While at low energies the only difference 
between the observed emission of the two jets will be as a 
result of the effects of beaming, at energies near the syn- 
chrotron cut-off the details of the acceleration mechanism 
will amplify this difference, depending on viewing angle. 

Although the work in this paper is limited to an ide- 
alised form of diffusion, and mildly relativistic shocks, it 
illustrates previous unexamined features which could be im- 
portant in the modelling of relativistic, 7-ray sources such 
as microquasars, blazars and GRBs. We have parameterised 
the exponential shape of the distribution cut-off and identi- 
fied new pitch angle dependent features between break and 
cut-off frequencies. Further work is needed to examine both 
momentum dependent scattering and high Lorentz factor 
flows. 



Figure 23. Synchrotron emission from the particle distribution 
shown in 1201 measured in the downstream medium. In plotting 
our fi = ±1 we used fj, = ±.9999 as there is no emission from an 
ordered field along fj, = ±1. 



Paniso, due to relativistic effects; and a cut-off, p cu t, due 
to energy losses. Given that the magnetic field is constant 
throughout this region it is trivial to produce the associ- 
ated synchrotron emission plots of figures 1211 1221 and 1231 It 
should be noted that in the emission plots I v is measured in 
the downstream frame, but since 3 is a Lorentz invariant 
the transformation is trivial. The synchrotron emission also 
includes the same three features we observed in the particle 
distribution; namely a break frequency beyond which the ef- 
fect of synchrotron cooling becomes important, a frequency 
at which pitch- angle or anisotropic effects play a role and an 
upper cut-off beyond which there is virtually no emission. 



4 DISCUSSION 

Particle acceleration and self-consistent synchrotron radia- 
tion have been considered previously by I Kirk et all (|l998l ) 
using a zonal model. They were successful in explaining the 
radio to X-ray spectrum of Mkn 501. However such zonal 
models typically depend on isotropic particle distributions. 
We have shown, however, that for particles near the high 
energy cut-off this is not true even for mildly relativistic 
flows. The computational resources available restricted our 
results to be below .7c. However even for the mildly rela- 
tivistic shock velocities we see a clear pattern of high energy 
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and here we present only Salzer summation: f(t, M) — > f(t) 
where 



This paper has been typeset from a TjrjX/ ]OTgX file prepared 
by the author. 



APPENDIX A: INVERSE LAPLACE 
TRANSFORMS 



While iHeavens fc Meisenheime'rl (|l987l) invert the Laplace 
transform analytically for particular cases here we use nu- 
merical methods as we will need to when dealing with rela- 
tivistic flows. 

Formally the inverse Laplace transform is the Bromwich 
integral, which is a complex integral given by: 



Kt) = c-^F {8)] = ±- i f 



■7 + ioo 



e 3t F(s) ds 



(Al) 



where 7 is to the right of every singularity of F(s). If the sin- 
gularity of F(s) all ly in the left half of the complex plane 
7 can be set to and this reduces to the inverse Fourier 
transform, which is easy to do. However for complicated or 
numerical Laplace functions the Bromwich integral is ex- 
tremely difficult to solve. The four main numerical inversion 
techniques are Fourier Series Expansion, Talbot's method, 
Weeks method and methods based on the Post-Widder for- 
mula. However some of these methods converge rather slowly 
and a lot of work has gone into creating acceleration meth- 
ods. Numerical Laplace inversion is a area of active research 
and the choice of inversion technique is as much an art as 
a science at the moment. In this paper the Post-Widder 
based method was chosen and only these methods shall by 
desc ribed below. Le t F(s) be the Laplace transform of f(t) 
then lWiddeil (] 19321 ) showed that f n (i) -> f(t) where 



/»(*) = 



(-ir 



^±ij /W((n + l)/i). (A2) 



The advantage of this method in our case in that we see 
that the Laplace transform of the solution times the Laplace 
coordinate is the zeroth order approximation to the actual 
solution. 



w> GO 



W(l/y)=>W {p)=pW(p). 



When dealing with numerical resu lts however it is easie r to 
use the Gaver-Stehfest algorithm (|Abate fc Varkdll2004l ). It 
is an algorithm based on the Post-Widder method with the 
Caver approximants, {f n (t) : n 0}, defined as 



/»(*) = 



(n + l)ln(2) (2(n + l) 



+ 1 



n+l 

£ 

fc=0 



n + l 
k 



/((n + l + fc)ln(2)/t). 



(A3) 



However the convergence for both these methods is slow. 
A test of methods for acc e lerati ng this convergence can be 
found in IValko fc Abate! (|2004l ) and two are found to be 
quite good: the non-linear Wynn's Rho Algorithm and the 
linear Salzer summation. Again a choice has to be made 



and 



f(t,M)=J2W k fr-i(t) 



W k = (-1)* 



M! \ k 



(A4) 



(A5) 



The Post-Widder method based on differentiation was im- 
plemented in Maple with the Salzer acceleration. It was used 
to produce the results in the non-relativistic limit as we 
have an analytic form of the Laplace function to work with. 
The Salzer accelerated Gaver-Stehfest algorithm was imple- 
mented in C/C++ code for use with the numerical output 
from the relativistic approach discussed above. 



APPENDIX B: DERIVING THE 
EIGENSYSTEM DIFFERENTIAL EQUATIONS 

The solutions, Qi, to equation 1361 

|- (^( M )(l- M 2 )^ -x\g( f j l )Q i = A i (u + »)Q i (Bl) 

for real x, are orthogon al, with we ight u + fj, and have real, 
distinct eigenvalues A; (Boas 1983). Taking the derivative of 
this equation with respect to x gives 



d 



D{n){l-n 2 ) 



d dQi 



dQi 



du \ cfyi dx J dx 

= A i ( u + ri^ + (^i(u + ») + \g((xfjQ i . (B2) 

Since solution to Sturm Liouville equations form an orthog- 
onal basis, we can write 



dQi _ ST n 



(B3) 



which gives 

^ q / d \ > 

E g " l 7T ( D ^( 1 ~ ' Lt2 )g _< 3 m ) - q m xXg(n)Qr, 

= E 9mA m (u + fl)Qr, 



= ^2 q m Ai(u + fj,)Q m + f -^-{u + a*) + ^g(fJ-) ) Qi 



dx 



(B4) 



Multiplying by Q* and integrating over a gives 
22 g m (A m — Ai){Qj, Qm) 

ill- 

= ^-(Qi>Qi) + x f gWQiQj (B5) 

Taking j = i implies equation [41] and j 7^ i implies equation 
HO] 



